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Abstract. The stability properties of inviscid protostellar disks are examined taking into account the Hall effect and buoyancy. 
Depending on the parameters, different types of instabilities can exist in different regions of disks. In a very low ionized region, 
the instability associated with baroclinic effects of buoyancy is likely most efficient. The Hall-driven shear instability can 
lead to destabilization of regions with a higher ionization. The magnetorotational instability modified by buoyancy can only 
be a destabilizing factor in regions with strong magnetic field or a relatively high conductivity (aB^ /p > Ueil, with ae the 
magnetization parameter of electrons). 
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1. Introduction 

The models of astrophysical disks require sufficiently strong 
turbulence to enhance the efficiency of angular momentum 
transport. At present, there is no commonly accepted view 
point as to how a laminar flow in disks is disrupted and turbu- 
lence generated. In general, turbulence may be generated due to 
various hydrodynamic and hydromagnetic instabilities which 
can arise in differentially rotating non-uniform gaseous disks 
but the exact origin of turbulence is still controversial. 

In accretion disks, the origin of turbulence is oftenly at- 
tributed to the magnetorotational instability since the necessary 
condition of instability, dil/dR < (i.e. a decrease of the 
angular velocity with cylindrical radius) is fulfilled (Velikov 
1959; Kurzweg 1963; Balbus & Hawley 1991; Kaisig, Tajima 
& Lovelace 1992; Kumar, Coleman & Kley 1994; Zhang, 
Diamond & Vishniac 1994). This instability exists not only 
for short wavelength perturbations, but also for global modes 
with scales comparable to the disk height (Curry, Pudritz & 
Sutherland 1994; CuiTy & Pudritz 1995). Note that the mag- 
netic shear instability can arise only if the field is not too 
strong, because this would suppress the instability (Urpin 
1996; Kitchatinov & Riidiger 1997). Simulations of the mag- 
netorotational instability in disks (Hawley, Gammie & Balbus 
1995; Matsumoto & Tajima 1995; Brandenburg et al. 1995; 
Torkelsson et al. 1996; Arlt & Riidiger 2001) show that the gen- 
erated turbulence may enhance the angular momentum trans- 
port. 

Unlikely, that the magnetic shear instability is the only in- 
stability that can exist in such complex objects as astrophysical 
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disks. A detailed analysis of MHD modes in stratified magnetic 
accretion disks demonstrates a much wider variety of instabil- 
ities than previously realized (Keppens, Casse & Goedbloed 
2002). Therefore, the current view on the origin of turbulence 
can be very simplified. Even a pure hydrodynamic origin of 
turbulence cannot be excluded (see Richard & Zahn 1999). 

The situation is particularly uncertain in cold and dense 
protostellar disks where the electrical conductivity is ex- 
tremely small because of a low ionization degree. The mag- 
netic Reynolds number is likely not very large in these disks, 
and the magnetic field cannot be considered as "frozen" into 
the gas (Gammie 1996). The behavior of the magnetic shear 
instability in the presence of ohmic dissipation has been con- 
sidered in the linear (Jin 1996) and non-linear regimes (Sano, 
Inutsuka & Miyama 1998). As it was first pointed out by 
Wardle (1999), however, poorly conducting protostellar disks 
can be strongly magnetized if electrons are the main charge car- 
riers. Magnetization leads to anisotropic electron transport with 
substantially different properties along and across the magnetic 
field (see, e.g., Spitzer 1978). If the field is sufficiently strong 
then the main contribution to the electric resistivity tensor is 
provided by the Hall component that produces the electric field 
perpendicular to both the magnetic field and electric current. 
This component is non-dissipative but it can change a geometry 
of the magnetic field. A linear stability analysis undertaken by 
Wardle (1999) shows that the Hall effect can provide an addi- 
tional either stabilizing or destabilizing influence depending on 
a direction of the magnetic field. A more general consideration 
of the magnetic shear instability in the presence of Hall cur- 
rents has been done by Balbus & Terquem (2001). They found 
that the Hall effect changes qualitatively the stability properties 
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of rotating gas resulting in destabilization of even outwardly 
increasing differential rotation (see also Riidiger & Shalybkov 
2003). In their analysis, however, Balbus & Terquem (2001) 
neglected gravity which, in protostellar disks, is as important 
as rotation. Gravity influences the behavior of modes via buoy- 
ancy, and may lead to new Hall-driven instabilities missed in 
simplified considerations. 

In the present paper we consider the linear stability prop- 
erties of magnetic protostellar disks taking into account the 
Hall effect and gravity. We treat the behavior of different short 
wavelength magnetohydrodynamic modes which can exist in 
such objects and determine the parameter domain where these 
modes are unstable. The paper is organized as follows. In 
Section 2, we discuss the Hall effect in the conditions of pro- 
tostellar clouds. In Section 2, the main equations are presented 
and a dispersion relation is derived that describes the behavior 
of short wavelength perturbations in the Boussinesq approxi- 
mation. The stability criteria for different modes are discussed 
in Section 4, and the growth rates of instabilities are calcu- 
lated in Section 5. Finally, our results are briefly summarized 
in Section 6. 

2. Anisotropic electric resistivity in protosteiiar 
disl<s 

The electrical conductivity is likely very low in protostellar 
disks because of a low temperature. The magnetic field can- 
not be considered as "frozen" into gas, and dissipative effects 
should be taken into account. However, despite low tempera- 
ture and ionization, the electron gas can be magnetized as it 
was first pointed out by Wardle (1999). Magnetization of the 
electron gas is characterized by the product of electron gy- 
rofrequency, lob = eB/nrieC and the relaxation time of elec- 
trons, T (see, e.g., Spitzer 1978). In protostellar disks, r is 
likely determined by the scattering of electrons on neutrals, 
then T = \/n{av) where {av) is the average product of the 
cross-section and velocity and n is the number density of neu- 
trals. Using the estimate of {av) obtained by Draine, Roberge 
& Dalgarno (1983) for electron-neutral collisions, we can rep- 
resent the magnetization parameter of electrons, a^, as 

Oe = UJBT = 21 Bn^lT-^'\ (1) 

where B is measured in Gauss, 7114 has units 10^^ cm^'', and 
T2 has units 100 K. If this parameter is greater than 1, i.e. 

B > 0.048 ni4 Gauss (2) 

(see Fig. then electron transport is anisotropic, and we 
have to use a tensorial magnetic diffusivity instead of a scalar 
one. For more details concerning the generalized Ohm's law 
in weakly ionized plasma we refer to the paper by Shalybkov 
& Urpin (1995) where this law has been considered using the 
relaxation time approximation for three components plasma 
with ions and neutrals of the same mass. In very weakly ion- 
ized plasma of protostellar disks, the difference between par- 
allel and perpendicular resistivity is small (see, e.g., Balbus 
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Fig. 1. The condition = 1 where the lines are marked with 
the values of the number density 77,14. For magnetic fields above 
the lines the Hall effect dominates the ohmic dissipation. 

& Terquem 2001). The Hall-originated magnetic diffusivity is 
then given by 

cB 

a.rj^pB^- , (3) 

An eric 

where ry is the microscopic magnetic diffusivity and rio is the 
number density of electrons. The magnetic diffusivity is 

rj = 2.34 X 10^ f^^T^^^ cm^ (4) 

where / = ric/n is an ionization fraction. 
The induction equation reads 

dB 

— = T0t(U xB) + riAB - /3rot(rotB x B). (5) 

The last term represents the Hall effect, and numerical evalua- 
tion done by Wardle & Ng (1999) indicate that this term can be 
of importance in some regions of protostellar disks. In (|5}, we 
neglect a non-uniformity of the resistivity and electron number 
density. 

3. Basic equations and thie dispersion relation 

Consider the stability properties of a magnetized axisymmet- 
ric protostellar disk of a finite vertical extent. The unperturbed 
angular velocity can generally depend on both R and z, so 
ri = ri(i?, z), where (R, cj), z) are cylindrical coordinates. The 
magnetic field, B = {Br, B^, B^), is assumed to be weak in 
the sense that the Alfven speed, ca, is small compared with 
the sound speed, Cg. This enables us to employ the Boussinesq 
approximation. 
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In the unperturbed state, the disk is assumed to be in hydro- 
static equilibrium in the R- and z-directions, 



Vp 



= G 



1 



YotB X B , G = g + n'^R , 



(6) 



where g is the gravity. If Cg > ca, the unperturbed Lorentz 
force is small compared with the pressure force, thus the disk 
structure is mainly determined by the balance between gravity, 
centrifugal force, and pressure. 

We consider axisymmetric short wavelength perturbations 
with the space-time dependence cxp(7t — ik ■ x) where k = 
(fcjjjOjfcz) is the wave vector, |fe • a;| ^ 1. Small perturba- 
tions will be indicated by subscript 1, whilst unperturbed quan- 
tities will have no subscript, except for indicating vector com- 
ponents when necessary. The linearized momentum and conti- 
nuity equations governing the behavior of such perturbations in 
the Boussinesq approximation read 

jU' + 2nxU' + e^R{U' ■ W)n = !^ _ aGT'+ 



— [k{B-B')-B'{k-B) 



k-U' ^0, 



(7) 
(8) 



where U', B' , p' and T' are perturbations of the hydrody- 
namic velocity, magnetic field, pressure and temperature, re- 
spectively; a — —{d\np/dT)p is the thermal expansion coef- 
ficient and 60 is the unit vector in the azimuthal direction. In 
Eq. (6) it is assumed that the density perturbation in the buoy- 
ancy force is determined by the temperature perturbation, thus 
p' = —paT', in accordance with the idea of the Boussinesq 
approximation. 

Since the thermal conductivity of protostellar clouds is low 
because of a low temperature (T ~ 10 — 10^ K), we adopt 
the adiabatic equation to describe the evolution of temperature 
perturbations. 



-fT' + U' ■ (AVT) = , 



(9) 



where (AVT) = VT - Vad?" is the difference between the 
actual and adiabatic temperature gradients. 

The linearized induction equation and the divergence free 
condition read 

(7 + ujk)B' = -iU'{k ■ B) + Re^{B' ■ Vf7) + 

+P{k-B)kxB', (10) 
k-B'^0, (11) 

where ljr = rjk'^ is the inverse timescale of the ohmic field 
decay. 

The dispersion equation governing the behavior of pertur- 
bations may be obtained in the standard way. Equating the de- 
terminant of the set of Eqs. 0. . . (^Jl to zero, we obtain 

7^ + 047^ + 037^ + a27^ + ai7 + ao = , (12) 

where 



a2 = 2LUn{Lul+u;l + Q'^), 

k 



and 

Q2 = 4f^2M 



2r!:^(2wH+Wsh) 



2nR 



dR 



ujI = -aAVr ■ 



• B) , 

R [Bn , dVL 



ojg is the frequency of buoyancy waves; = {k ■ B) j \J\-Kp is 
the Alfven frequency, represents the effects associated with 
differential rotation, wh and are the characteristic frequen- 
cies of the Hall- and shear-driven processes, respectively. 

If gravity is neglected (ojg = 0) then we recover the disper- 
sion equation derived by Balbus & Terquem (2002) 
3 



where 

63 = 2wr, 



0, 



(13) 



+ wh(wh + t^sh) + 2^1 + ' 



&1 ^ 2tjR(tji + g2 



fen = < 



k 



Assuming that the condition of instability is given by 

&o <0, 



(14) 



Balbus & Terquem (2001) argued that the Hall effect can desta- 
bilize any differential rotation laws in protostellar disks, even 
those with angular velocity increasing outward. 

4. Criteria of instabiiity of protosteiiar disi<s 

The equation ( I12> describes five low-frequency modes which 
can exist in protostellar disks. The condition that at least one 
of the roots of equation J12t has a positive real part (unstable 
mode) is equivalent to one of the following inequalities 

ao < , 

A\ = 0403 — 02 < , 

A2 = 02(0403 — 02) - 04(0401 — oo) < , 

Az = (04O1 - 00)[02(0403 - 02) - 04(0401 - oo)]- 



-00(0403 - 02)^ < , 



(15) 



being fulfilled (see, e.g., Aleksandrov, Kolmogorov & 
Laurentiev 1985). Since cjr is positive defined quantity, the 
first condition (04 < 0) will never apply, and only the four 
given conditions determine the instability in disks. 
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4. 1. The condition oq < 

Since u!\ > 0, the condition ao < is equivalent to 

ujI < 0, 



or 



fc2(G • AVT) - (fc • G)(fe • AVT) > 0. 



Generally, tjg < if the temperature gradient is superadiabatic, 
i.e. it exceeds its adiabatic value. In this case, the standard con- 
vective instability arises. However, may also be negative if 
the temperature gradient is subadiabatic but AVT is not par- 
allel to the "effective gravity", G (see Urpin & Brandenburg 
1998). This obliqueness can be caused, in principle, either by 
the dependence of fi on z or by radiative heat transport in the 
radial direction. Introducing the angle ip between the vectors G 
and k and representing AVT as a sum of components parallel 
and perpendicular to G, AVT = (AVr)|| + (AVT)^, the 
inequality (I17> can be rewritten in the form 

G(AVr)||[sin2 V'-sinV'cosV'(AVr)_L/(AVT)||] > .(18) 

If stratification is stable according to the standard 
Schwarzschild criterion of convection, G(AVT)|| < 0, 
then the instability arises at 



sin2V-sinV'C0sV'(AVT)i/(AVr)|| <0. 



(19) 



This condition can be fulfilled due to the obliqueness of G and 
AVT for perturbations with a small (but non-zero) angle ip. 
Estimating (AVT)^ as (z/i?)(AVT)||, we obtain that w| for 
unstable perturbations is small, 



(20) 



Therefore, convection caused by obliqueness of G and AVT 
is probably relatively slow. 

4.2. The condition Ai <Q 

This condition reads 



'^A + '^R + + t^H^sh < 0, 



or, substituting the frequencies. 



(3R{k-B) [k. 



dR 



dQ 
1): 



(21) 



(3^k\k-BY. (22) 



Despite this inequality is the only condition (I15> that does not 
depend on gravity, it differs from the criterion of instability de- 
rived by Balbus & Terquem (2001) (see Eq. (85) of their paper). 
The condition M2\ describes a new instability that appears due 
to combined influence of shear and the Hall effect. This insta- 
bility differs from the magnetic shear instability because the 
only term that can provide a destabilizing influence is propor- 
tional to the Hall frequency and shear stresses, and this terms 
is vanishing if lo-h 0. To satisfy the inequality M2\ the sign 
of the left hand side should be negative since all three terms on 
the right hand side are negative. Obviously, for any dependence 



of O on i? and z and for any direction of B, there exist wave 
vectors that satisfies the inequality 



(16) {k-B)[ k 



(17) 



OR 



kR 



dz 



< 0, 



(23) 



and makes the left hand side of Eq. i22\ negative. Therefore, 
any differential rotation can generally be unstable if the Hall 
effect is sufficiently strong. 

Compare characteristic frequencies in Eq. (I2H . If the Hall 
parameter is large, Ce S> 1, then the "ohmic frequency" is 
negligible in Eq. i21\ . The characteristic value of the Hall fre- 
quency, LUn, is 



(24) 



where A = 27r/fc is the wavelength, and An has units 10^^ cm; 

are in units 100 cm^'^. Assuming that gas is weakly ionized 
and n ^ rio, we can estimate 



LUA 



ck 



1.4 X io-4Ar>;//^/ 



1/2 r-1/2 



(25) 



where ujp — {ine^nc/mcY^'^ is the plasma frequency, rric 
and TOp are the electron and proton mass, respectively. Since 
/ can be very small in a large fraction of the disk volume 
(/ 10"("- -i3)), most likely that |wh| > |wa|- In this case, 
the Alfven frequency is small in Eq. ( I2H . Then, in a row with 
Eq. i23i . the necessary condition of instability is also 



(26) 



|Wsh| > |wh|. 



Assuming \dn/dR\ ^ Vl/R as is usual in astrophysical disks, 
we can represent the necessary condition \26\ as 



B < 



10 -^P /nc2Afi Gauss, 



(27) 



where P = 27r/i7 is the rotation period, Py^ has units 1 yr. 
The condition that the Hall effect dominates ohmic dissipation, 
fle > 1, yields 



B > 0.048 nu T^^^ Gauss. 

The conditions ( I27> and ( I28> are consistent only if 



Py, < 2.2 X 10-2/_i2A2iT2 



(28) 



(29) 



where /_i2 = //lO^^^. If we estimate the thickness of a disk, 
H, at the distance i? as 7? ~ O.li? and assume that A ^ O.IH 
to justify a short wavelength approximation, then we have at 
i? = 1 AU, 



Py, < 2.2 X 10-2/_i2T2 



(30) 



To be fulfilled this condition requires the temperature lower 
than 100 K and ionization higher than 10^^^. Therefore, the 
condition i29\ seems to be very restricting and hardly to be 
fulfilled in the conditions of protostellar disks. 
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4.3. The condition ^2 < 

Substituting the expressions for coefficients, we can represent 
the condition A2 < as 



(31) 



All terms on the left hand side are positive except and, 
hence, the criterion i31\ can be satisfied only if < 0. 
However, this criterion requires larger negative lo^ than the 
condition J16t . Therefore, the criterion ^2 < can be fulfilled 
only if the condition ao < is already fulfilled. 



4.4. Ttie condition A3 < 

This condition generalizes the criterion obtained by Balbus and 
Terquem (2001) for the protostellar disk with no gravity. If we 
assume g — (and, hence, uj^ = qq = 0) then A3 — a4,aiA2- 
Since 04 is always positive and A2 is positive at g = (see Eq. 
J31» . the criterion A3 < reduces to 



ai|g=o = &o < 0, 



(32) 



that is the condition derived by Balbus & Terquem (2001). 
In the general case, the criterion A3 < yields 



D = bo + q{[u 



+0)1 + 2 ( 2f7^ 



<0, 



where 
1 



(21]% 



(33) 



(34) 



Gravity can provide either positive or negative contribution to 
the left hand side of Eq. ( I33> and, hence, can be either a stabiliz- 
ing or destabilizing factor depending on the characteristic fre- 
quencies. The condition (I33> is generally rather cumbersome, 
and we consider only the particular cases of astrophysical in- 
terest. 

In protostellar disks, we have typically Wg ~ and, hence, 
q ^ 1. As it has been adopted in previous studies, we also 
assume that the Hall effect dominates ohmic dissipation, and 

Wh > tJR (or fle > 1), 



c^R « 9.2 X 10-Vri2T2^'Ari' s 



-2 „-i 



(35) 



Stability properties are sensitive to the relationship between the 
Hall frequency, cjh, and the angular velocity, il. The inequality 
Wh > is fulfilled only if the field is strong enough. 



(36) 



This field seems to be rather strong for the conditions of pro- 
tostellar disks and, most likely, ^ o^h. Note also that at 
wh > we have D > for both cjh > and cja > t^H- 



Hence, the instability does not appear, and the field, satisfy- 
ing Eq. ( I36> . is stabilizing for instability given by the criterion 
A3 < 0. Therefore, we consider Eq. ( I33t only at ^ wh- 
If this is the case, the criterion of instability reads 

t^H(t^H + t^sh) [(1 + g)(^g + Q^) - + + 



2r!y(c^H 



+ 



2"« 



0' < 0. 



(37) 



As it was mentioned (see Eq. ( I25» . the Hall frequency can ex- 
ceed the Alfven frequency in protostellar disks, cjh > t^A, and 
we treat this case first. Under this assumption, the criterion ( I37> 
yields 



(38) 



In a convectively stable disk with uj^ > and > 0, the 
shear-driven instability arises only if 



(39) 



To be fulfilled, this inequality requires a strong shear, satisfy- 
ing Eq. (I26> . and a wave vector directed in accordance with Eq. 
(I23> . There is no instability due to the Hall effect without shear. 
It appears that the condition ujsh > i^H is general for instabili- 
ties caused by a joint influence of the Hall effect and shear but, 
as we akeady mentioned, it is very restricting in protostellar 
disks. The sign of ujsh plays no crucial role in the condition 
( I39> , and instability can appear for any dil/dR if Eq. (I26> is 
fulfilled. Note that in the case g = 0, the instability can arise 
for any sign of shear as well (Balbus & Terquem 2001). 

Consider now the criterion of instability (I37> in the case 
when the Hall frequency is small compared to the Alfven 
frequency, uja > ^h- We again assume il > cjh but the 
Alfven frequency can generally be larger or smaller than O. 
If U!A ^ or 



B> 0.15 nli^ An P- 



yr ) 



(40) 



then D sa u!\ > 0, and instability is not possible. The field, sat- 
isfying this condition, is typically even stronger than that given 
by (I36> . Unlikely that such a strong field can exist in proto- 
stellar disks, and we consider instability under a more realistic 
condition ft ^ o^a- 

If n > UJA > ^ yc^n/n, then we have from the condition 

(1 + 9) 



< 0. 



In a convectively stable disk, the instability appears if 



(41) 



(42) 



At cjg ~ 51 ^ UJA, this condition is far more restricting than 
the criterion of the magnetorotational instability obtained by 
Balbus & Terquem (2001) for protostellar disks. 



(El}. 



If ily^oju/il > UJA, then we recover the instability criterion 
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5. The growth rate of instability 

A general analysis of the roots of Eq. lll2t is very complicated. 
However, simple expressions for the roots can be obtained in 
some cases of astrophysical interest. Consider initially Eq. il2\ 
in the case when the ohmic dissipative "frequency", ujn, is 
small compared to other characteristic frequencies. Then, with 
accuracy in terms of the zeroth order in ojr, four roots cor- 
responding to rapidly varying modes (either growing, or de- 
caying, or oscillating) can be obtained from Eq. il2\ with ne- 
glected ohmic dissipation. These roots satisfy the approximate 
equation 



7'' + C27^ + Co = 0, 
where 



(43) 



(44) 



Co = ioi [uj'g + 2n—ujshj + l^Mloh + c^siOK + Q )■ ^"^5) 

The fifth mode describes a secular instability and varies on a 
long time scale proportional to lu^. For this mode, we have 
with the accuracy in terms of the lowest order in wr. 



75 



Oq 

ai 



1 

Co 



(46) 



If uj\ 3> IwHC'jJH+t^sh)!, then Eq. (46) yields a simple estimate 



for 75, 



44 fzloT^'^XT?P. 



12^2 ^ni yr- 
The solution of Eq. ( I43t is given by 



7 



C2 



± 



Co- 



(47) 



(48) 



Note that the true expressions for 7 have to contain also small 
corrections, A7, proportional to cjr which are neglected in Eq. 
(48). If 17 > max(wH, wa) as the criteria of instability requires 
then C2 ^ 4co, and we obtain 



1I2 



Co 
C2' 



73% 



-C2- 



(49) 



The modes 1 and 2 are relevant to the Hall-driven instability 
whereas the modes 3 and 4 describe convection and are stable 
in convectively stable disks. One of the modes 71 2 is unstable 
if Co < 0. Note that the secular mode ( I46> is unstable as well 
under this condition, and there exist two qualitatively different 
instabilities represented by the same criterion. 

In the case fi-^/wn/i^ > ^^A, we have for the Hall-driven 
modes 



7l 2 ~ -t^H(wH + Wsh)- 



(50) 



If the condition ( I39> is fulfilled then one of these modes is un- 
stable. Since Wsh ~ we can estimate 



71,2 



^ 31 SV2„-V2 



^11 -'^yr • 



(51) 



If > OJA > ^ v^j^h/^ then we obtain 



I1.2 



(52) 



This is the known dispersion equation for magnetic shear- 
driven modes (see, e.g., Urpin 1996) in the limit ujg ^ oja, 
and one of these modes is unstable if 



2n—ujsh < 



(53) 



The growth rate is of the order of uja- Note that the condition 
(53) differs from the more general condition ( I42> that predicts 
instability also within the range 



4 < 2n^uj,^ < -^uil 



(54) 



fc"°" ~ 2 

The difference originates from the approximate expressions 
( I49> that describes 7J 2 neglecting terms proportional wr. If 
the inequality ( I53t is not fulfilled, we have cq > 0, and the 
roots 71 and 72 become imaginary. Therefore, their stability at 
Cq > is determined by small real corrections cx wr. These 
corrections can be calculated by making use of the standard 
perturbation procedure. For the modes 1 and 2, corrections are 
same and given by 

A7i,2 - -^(ao-2a4Co) = ( ^n'^uj^h + c^l) .(55) 



2co 



These corrections are positive under the condition i54\ and 
leads to instability in complete agreement with the criterion 
(I42> . However, the instability caused by small resistive correc- 
tions is oscillatory and qualitatively different from the magnetic 
shear instability. Note that within the range ( I54> the secular 
mode i46i is stable. 

Consider now the roots of Eq. (I12> assuming that the dis- 
sipative "frequency", wr, is large compared to other charac- 
teristic frequencies including cjh- Most likely, that the angular 
velocity, il, is largest among other frequencies, and the assump- 
tion wr > fl implies that 

1/2 



Py,. > 2.2 X 10 



'■/-12A11T2 



(56) 



Probably, such inequality can be fulfilled in some regions of 
protostellar disks if ionization is low. If > oja, then the co- 
efficients of Eq. J12t are given by 



04 w 2tjR, 03 w Wr, 02 w a4(Wg 
ai w ujRa2/2, ao « WRwi^g. 



Four roots of the Eq. (I12t with such coefficients have negative 
real parts, and only one root can correspond to instability. With 
the accuracy in terms of the lowest order in ijJr^, this root is 
given by 



7 



wr uj. 



(57) 



The growth rate is positive if ut"^ < 0. As it was argued 
above (see Eq. il6\ ). the instability can appear under this con- 
dition if G and AVT are not parallel. In protostellar disks, 
the growth rate of this instability is probably relative small, 
7 ^ ll!a{'^aI^'r){H / Kf' . This instability, however, can ex- 
ist if ionization is extremely low, and the magnetic Reynolds 
number is small even for rotation. 
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6. Conclusions 

We have considered the stability properties of protostellar 
disks. Compared to other studies we took into account con- 
sistently gravity and buoyancy forces which play an important 
role in dynamics of disks. It has been shown that, under certain 
conditions, magnetic protostellar disks are unstable to different 
kinds of shear-driven instabilities. Some of these instabilities 
are relevant to the Hall effect but some can manifest themselves 
even if the Hall currents are negligible. All considered instabili- 
ties are influenced by the magnetic field, and a relatively strong 
field can suppress instabilities. The strength of the field which 
can stabilize the flow depends generally on the conditions in 
disks and on the wavelength of the perturbation. The stabiliz- 
ing field is determined by the conditions i36\ and ( I40> and is 
typically relatively strong. 

The type of instability that is more efficient depends very 
much on the conditions in protostellar disks and the wavelength 
of perturbations, and can be different in different regions of the 
disk. The instability associated with the criterion ( I16t seems to 
be most general in protostellar disks. It can occur both in mag- 
netized (oe > 1) and non-magnetized (oe < 1) regions and 
requires only non-parallel G and AVT. The growth rate of 
this instability is given by equations (46) and (57) in the lim- 
iting cases of small and large ujj^, respectively. Note that the 
baroclinic instability appears if dissipative processes are taken 
into account and, therefore, it can also exist in accretion disks. 
In the present study, the baroclinic convection is caused by the 
electrical resistivity and magnetic field but it can also exist, for 
example, in non-magnetic disks if one takes account of the ther- 
mal conductivity (see Urpin & Brandenburg 1998). The baro- 
clinic convection, however, can be not the most efficient insta- 
bility operating in protostellar disks since the inequality il6\ is 
fulfilled only for a very narrow range of k almost parallel to 
G. In magnetized regions with > 1, other instabilities can 
manifest themselves as well if the wavelength of perturbations 
is sufficiently long. For the purpose of illustration, we consider 
the domains of different instabilities in the case when the wave 
vector is not parallel or perpendicular to Jl and B and, hence, 
kz ^ k and {k ■ B) ^ kB. Then, the regions of instability 
are determined by two characteristic wavelength. Ah and Aa, 
corresponding to the conditions o^h = ^ and loj^ = fl, respec- 
tively. These wavelengths are 

Ah = 3.1 X 10^^Si/2„-^i/2pV2 (58) 

Aa = 6.6 X 10"Sn~4 /^Pyr cm, (59) 
and their ratio is given by 

Ah 

If ^ ^ 1 then for A ^ Aa (that is equivalent to 51 ^ uja ^ 
ft y/iJii/Tt) a particular case of the magnetic shear instability 
may occur with the growth rate given by Eq. (52). If f <C 1 and 
A ^ Ah (that is equivalent to £7 ^ fl o^h / ft ^ wa) then the 
Hall-driven shear instability can arise with the growth rate (50). 
Note that the both these instabilities co-exist with the baroclinic 
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Fig. 2. The same as in Fig.^but only for = 1 

convection in their domains of instability but the latter can also 
occur for shorter A. 

In the present paper, we have addressed the behavior of 
only axisymmetric perturbations. It is clear, however, that the 
obtained results can apply to nonaxisymmetric perturbations 
with the azimuthal wavelength much longer than the vertical 
or radial ones, min(fcfl, k^) ^ fc^. The turbulence that could 
be generated by the considered instabilities may be strongly 
anisotropic in the {R, z)-plane, because the instability criteria 
are very sensitive to the direction of the wave vector. However, 
the generated turbulence may be efficient in the radial transport 
of angular momentum. 
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